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We consider the effect of a small-scale matter-antimatter domain structure on big bang nucleosyn- 
thesis and place upper limits on the amount of antimatter in the early universe. For small domains, 
which annihilate before nucleosynthesis, this limit comes from underproduction of 4 He. For larger 
domains, the limit comes from 3 He overproduction. Most of the 3 He from p 4 He annihilation is anni- 
hilated also. The main source of 3 He is photodisintegration of 4 He by the electromagnetic cascades 
initiated by the annihilation. 
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If the early universe was homogeneous, antimatter an- 
nihilated during the first millisecond. However, baryoge- 
nesis could have been inhomogeneous, possibly resulting 
in a negative net baryon number density in some regions 
After local annihilation these regions would have 
only antimatter left, resulting in a matter- antimatter do- 
main structure. 

There are many proposed mechanisms for baryogene- 
sis [Q. In models connected with inflation, there is no 
a priori constraint on the distance scale of the matter- 
antimatter domain structure that may be generated. If 
the distance scale is small, the antimatter domains would 
have annihilated in the early universe, and the presence 
of matter today indicates that originally there was less 
antimatter than matter. 

We consider here such a scenario: a baryoasymmet- 
ric universe where the early universe contains a small 
amount of antimatter in the form of antimatter domains 
surrounded by matter. We are interested in the effect of 
antimatter on big bang nucleosynthesis (BBN) J3|. Much 
of the earlier work on antimatter and BBN [0-CO| has fo- 
cused on a baryon-antibaryon symmetric cosmology m 
or on homogeneous injection of antimatter through some 
decay process ||. 

The smaller the size of the antimatter domains, the 
earlier they annihilate. Domains smaller than 100 m at 
1 MeV, corresponding to 2 x 10~ 5 pc today, would an- 
nihilate well before nucleosynthesis and would leave no 
observable remnant. 

The energy released in annihilation thermalizes with 
the ambient plasma and the background radiation, if the 
energy release occurs at T > 1 keV. If the annihilation oc- 
curs later, Compton scattering between heated electrons 
and the background photons transfers energy to the mi- 
crowave background, but is not able to fully thermalize 
this energy. The lack of observed distortion in the cosmic 
microwave background (CMB) spectrum constrains the 
energy release occuring after T = 1 keV to below 6 x 10~ 5 
of the CMB energy pi |12|. This leads to progressively 



stronger constraints on the amount of antimatter anni- 
hilating at later times, as the ratio of matter and CMB 



energy density is getting larger. Above T ~ 0.1 eV the 
baryonic matter energy density is smaller than the CMB 
energy density, so the limits on the antimatter fraction 
annihilating then are weaker than 6 x 10~ 5 . 

For scales larger than 10 pc (or 10 11 m at T = 1 
keV) the tightest constraints on the amount of antimat- 
ter come from the CMB spectral distortion, and for even 
larger scales from the cosmic diffuse gamma spectrum 
& 

We consider here intermediate domain sizes, where 
most of the annihilation occurs shortly before, during, or 
shortly after nucleosynthesis, at temperatures between 1 
MeV and 10 eV. The strongest constraints on the amount 
of antimatter at these distance scales will come from BBN 
affected by the annihilation process. 

Rehm and Jedamzik J^] considered annihilation imme- 
diately before nucleosynthesis, at temperatures T = 80 
keV - 1 MeV. Because of the much faster diffusion of neu- 
trons and antineutrons (as compared to protons and an- 
tiprotons) the annihilation reduces the net neutron num- 
ber Q, leading to underproduction of 4 He. This sets a 
limit R < few % to the amount of antimatter relative to 
matter in domains of size ta ~ 1 cm at T = 100 GeV 
(4 x 10 6 m at T = 1 keV). 

We extend these results to larger domain sizes, for 
which annihilation occurs during or after nucleosynthe- 
sis. Since our results for the small domains and early 
annihilation agree with Rehm and Jedamzik, we concen- 
trate on the larger domains and later annihilation in the 
following discussion. Below, all distance scales given in 
meters will refer to comoving distance at T = 1 keV. 

The case where annihilation occurs after nucleosynthe- 
sis was considered in Refs. 0. Because annihilation of 
antiprotons on helium would produce D and 3 He they es- 
timated that the observed abundances of these isotopes 
place an upper limit R < 10 -3 to the amount of anti- 
matter annihilated after nucleosynthesis. As we explain 
below, the situation is actually more complicated. 

Consider the evolution of an antimatter domain (of di- 
ameter 2r) surrounded by a larger region of matter. At 
first matter and antimatter are in the form of nuclcons 
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and antinucleons, after nucleosynthesis in the form of ions 
and anti-ions. Matter and antimatter will get mixed by 
diffusion and annihilated at the domain boundary. Thus 
there will be a narrow annihilation zone, with lower den- 
sity, separating the matter and antimatter domains. At 
lower temperatures (T < 30 keV) the pressure gradient 
pij drives matter and antimatter towards the annihila- 
tion zone. This flow is resisted by Thomson drag, which 
leads to diffusive flow |L5|]. 

Before nucleosynthesis, the mixing of matter and an- 
timatter is due to (anti)neutron diffusion. When He 
is formed, free neutrons disappear, and the annihilation 
practically ceases. If annihilation is not complete by 
then, it is delayed significantly because ion diffusion is 
much slower than neutron diffusion. There will then be 
a second burst of annihilation well after nucleosynthesis, 
at T ~ I keV or below. Indeed, depending on the size of 
the antimatter domains, most of the annihilation occurs 
either at T > 80 keV (for r < 2 x I0 7 m) or at T < 3 
keV (for r > 2 x I0 7 m). 

The annihilation is so rapid that the outcome is not 
sensitive to the annihilation cross sections. The ex- 
act yields of the annihilation reactions are more impor- 
tant. From the Low-Energy Antiproton Ring (LEAR) at 
CERN, we have data for antiprotons on helium, and also 
for some other reactions with antiprotons p~8|] . 

The annihilation of a nucleon and an antinucleon pro- 
duces a number of pions, on average 5 with 3 of them 
charged jlTj]. The charged pions decay into muons and 
neutrinos, the muons into electrons and neutrinos. The 
neutral pions decay into two photons. About half of the 
annihilation energy, 1880 MeV, is carried away by the 
neutrinos, one third by the photons, and one sixth by 
electrons and positrons jl. 

If the annihilation occurs in a nucleus, some of the 
pions may knock out other nucleons. Part of the anni- 
hilation energy will go into the kinetic energy of these 
particles and the recoil energy of the residual nucleus. 
Experimental data on the energy spectra of these emitted 
nucleons are well approximated by the formula Ce~ E / E ° , 
with average energy Eq ~ 70 MeV, corresponding to a 
momentum of 350 MeV/c @. 

After 4 He synthesis, the most important annihilation 
reactions are pp and p 4 He. According to Balestra et al. 
fyl, a p 4 He annihilation leaves behind a 3 H nucleus in 
43.7±3.2 % and a 3 He nucleus in 21.0±0.9 % of the cases. 
The rms momentum of the residual 3 He was found to be 
198 ± 9 MeV/c. 

It is important to consider how these annihilation 
products are slowed down. If they escape far from the 
antimatter domain, they will survive; but if they are ther- 
malized close to it they will soon be sucked into the an- 
nihilation zone 0. 

Fast ions lose energy by Coulomb scattering on elec- 
trons and ions. If the velocity of the ion is greater than 
thermal electron velocities, the energy loss is mainly due 



to electrons. At lower energies the scattering on ions 
becomes more important. Below T — 30 keV, when 
the thermal electron-positron pairs have disappeared, the 
penetration distance of an ion of initial energy E depends 
on the ratio E/T [pL For E > (M ion /m e )T, the pene- 
tration distance is [M 
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where A ~ 15 is the Coulomb logarithm, giving a comov- 
ing distance 



1 1 (E 



comoving 



AZ 2 mo \T 



0.4 m. 



(2) 



For smaller E/T, I keeps getting shorter, but not as fast 
as Eqs. (@) and @ would give @. 

For 3 H and especially for 3 He, I would become compa- 
rable to the original size of the antimatter domain only 
well after the annihilation is over. Thus only a small frac- 
tion of these annihilation products escape annihilation. 
For D this fraction is larger, but still small, except for 
the largest domains considered here. 

Neutrons scatter on ions, losing a substantial part of 
their energy in each collision. The neutrons from anni- 
hilation reactions have sufficient energy to disintegrate a 
4 He nucleus. This hadrodestruction [5l| of 4 He causes 
some additional 3 He and D production. Because pro- 
tons are more abundant than 4 He nuclei, a neutron is 
more likely to scatter on a proton. The mean free path 
A = l/((T np n p ) is larger than the distance scales con- 
sidered here, so the annihilation neutrons are spread out 
evenly. At lower temperatures (T < 1 keV), neutrons de- 
cay into protons before thermalizing. At higher temper- 
atures, the stopped neutrons form deuterium with pro- 
tons. 

The high-energy photons and electrons from pion de- 
cay initiate electromagnetic cascades 21-23]. Below 
T = 30 keV, the dominant processes are photon-photon 
pair production and inverse Compton scattering 
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with the background photons ^b- The cascade photon 
energies E 1 fall rapidly until they are below the threshold 
for pair production, i? 7 e 7 = m 2 , where e 7 is the energy of 
the background photon. Because of the large number of 
background photons, a significant number of them have 
energies 3> T, and the photon-photon pair production is 
the dominant energy loss mechanism for cascade photons 
down to |E3| 
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When the energy of a 7 falls below E max , its mean free 
path increases and it is more likely to encounter an ion. 



2 



As the background temperature falls this threshold en- 
ergy rises, and below T ~ 5 keV, £ max becomes larger 
than nuclear binding energies, and photodisintegration 
pl| p4| becomes important. Photodisintegration of D 
begins when the temperature falls below 5.3 keV, photo- 
disintegration of 3 He ( 3 H) below 2.2 keV (1.9 keV) and 
photodisintegration of 4 He below 0.6 keV. 

Thus there are two regimes for photodisintegration: (1) 
between T = 5.3 keV and T = 0.6 keV, where the main 
effect is photodisintegration of D, 3 H, and 3 He; and (2) 
below T = 0.6 keV, where the main effect is production 
of these lighter isotopes from 4 He photodisintegration. 
Because of the much larger abundance of 4 He, even a 
small amount of annihilation during the second regime 
swamps the effects of the first regime, and only in the 
case that annihilation is already over by T = 0.6 keV, 
is D photodisintegration important [pl"|-|2~i|l . Because of 
the difference in the neutron and proton diffusion rates, 
domains this small have already significant (neutron) an- 
nihilation before 4 He synthesis. 

For the larger domain sizes, the most significant ef- 
fect of antimatter domains on BBN turns out to be 3 He 
production from He photodisintegration. 

We have done numerical computations of nucleosyn- 
thesis with antimatter domains. Our inhomogeneous nu- 
cleosynthesis code includes nuclear reactions, diffusion, 
hydrodynamic expansion, annihilation, spreading of an- 
nihilation products, photodisintegration of 4 He and dis- 
integration by fast neutrons pj| |. Because of the lack 
of data on the yields of annihilation reactions between 
nuclei and antinuclei, we have not incorporated antinu- 
cleosynthesis in our code, but the antimatter is allowed 
to remain as antinucleons. This could affect our results 
at the 10% level. 

For photodisi nteg ration of 4 He we use the results of 
Protheroe et al. J24J scaled by the actual local 4 He abun- 
dance. The 3 He yield is an order of magnitude greater 
than the D yield [gij . The Protheroe et al. results assume 
a standard cascade spectrum. This will not be valid for 
temperatures below 100 eV, since E m&K becomes compa- 
rable or greater than typical energies of the initial 7's 
from annihilation. It would be important to find out the 
true cascade spectrum for these low temperatures, since 
this will affect our results at the largest scales. 

We have in mind a situation where antimatter domains 
of typical diameter 2r are separated by an average dis- 
tance 2L. This we represent with a spherically symmet- 
ric grid of radius L with antimatter at the center with 
radius r. For simplicity we assume equal homogeneous 
initial densities for both the matter and antimatter re- 
gions. This density is set so that the final average density 
after annihilation will correspond to a given baryon-to- 
photon ratio 77. Since we are looking for upper limits 
to the amount of antimatter which come from a lower 
limit to the He abundance and an upper limit to the 
3 He abundance, we choose 77 = 6 x 10~ 10 , near the up- 



per end of the acceptable range in standard BBN, giving 
high 4 He and low 3 He. 

We show the 3 He yield as a function of the antimatter 
domain radius r and the antimatter/matter ratio R in 
Fig. |l|. For domains smaller than r = 10 5 m, annihila- 
tion happens before weak freeze-out, and has no effect on 
BBN. For domains between r = 10 5 m and r — 10 7 m, 
neutron annihilation before 4 He formation leads to a re- 
duction in 4 He and 3 He yields. For domains larger than 
r = 2 X 10 m, most of the annihilation happens after 
4 He synthesis. Antiproton-helium annihilation then pro- 
duces 3 He and D, but most of this is deposited close to 
the annihilation zone and is soon annihilated. The much 
more important effect is the photodisintegration of He 
by the cascade photons, since it takes place everywhere in 
the matter region and thus the photodisintegration prod- 
ucts survive. This leads to a large final 3 He and D yield. 
The same applies to n 4 He reactions by fast neutrons from 
p 4 He annihilation, but the effect is much smaller, because 
p 4 He annihilation is less frequent than pp annihilation, 
and a smaller part of the annihilation energy goes into 
neutrons than in the electromagnetic cascades. 
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FIG. 1. The 3 He yield as a function of the antimat- 
ter/matter ratio R and the antimatter domain radius r. The 
distance scales are given both at T = 1 keV (in meters) and 
today (in parsecs). We plot contours of (the logarithm of) 
the number ratio 3 He/H. The dotted lines show contours of 
the "median annihilation temperature", i.e., the temperature 
of the universe when 50% of the antimatter has annihilated. 
Typically the annihilation is complete at a temperature lower 
than this by about a factor of 3. 

We obtain upper limits to the amount of antimatter in 
the early universe by requiring that the primordial 4 He 
abundance Y p must not be lower than Y p = 0.22, and 
that the primordial 3 He abundance must not be higher 
than 3 He/H = 10~ 4 - 5 ||. (The standard BBN results for 
77 = 6 X 10~ 10 are Y p = 0.248 and 3 He/H = 1.1 x 10~ 5 .) 
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For domain sizes r < 10 11 m (or 10 parsecs today), these 
limits are stronger than those from the CMB spectrum 
distortion. See Fig. 0. 
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FIG. 2. Upper limits from BBN (solid lines) and CMB 
(dot-dashed line) to the antimatter/matter ratio J? as a func- 
tion of the antimatter domain radius r. The dashed line gives 
the upper limit from BBN if photodisintegration is ignored. 



For r < 10 5 m, there is no BBN constraint on antimat- 
ter. For r = 10 5 -10 7 m, the amount of antimatter can be 
at most a few per cent, to avoid 4 He underproduction. 
Our limit is somewhat weaker than that of Rehm and 
Jedamzik, since they considered a lower r\ = 3.4 x 10 -10 . 

For larger domains, antimatter annihilation causes 3 He 
production from He photodisintegration and the limit 
reaches R = 2 x 1CT 4 at r ~ 10 9 m. 

There may exist small regions of parameter space 
where acceptable light element yields would be obtained 
for "nonstandard" values of ij and large R Clearly 
the simultaneous reduction of 4 He and increase of 3 He 
and D suggest such a possibility for large rj. 
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